Nonsequential Double Recombination in Intense Laser Fields 
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A second plateau in the harmonic spectra of laser-driven two-electron atoms is observed both 
in the numerical solution of a low-dimensional model helium atom and using an extended strong 
field approximation. It is shown that the harmonics well beyond the usual cut-off are due to the 
simultaneous recombination of the two electrons, which were emitted during different, previous 
half-cycles. The new cut-off is explained in terms of classical trajectories. Classical predictions 
and the time-frequency analysis of the ab initio quantum results are in excellent agreement. The 
mechanism corresponds to the inverse single photon double ionization process in the presence of a 
(low frequency) laser field. 

PACS numbers: 42.65.Ky, 32.80.Rm, 34.80.Lx, 32.80.Wr 
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High-order harmonic generation (HOHG) is one of 
the fundamental processes that occur when intense laser 
pulses interact with atoms or molecules (see, e.g., Refs. 

Q for recent reviews) . In the case of a linearly polar- 
ized incoming laser field, odd multiples of the laser fre- 
quency are emitted with a relatively high and almost con- 
stant efficiency (~ 10~^) up to the celebrated /p-|-3.17 Up 
cut-off 3] , where /p is the ionization potential and U-p is 
the ponderomotive energy (i.e., the time-averaged quiver 
energy of a free electron in the laser field). In such a way, 
coherent short-wavelength radiation down to the "water 
window" could be generated using "table-top" equipment 

HOHG up to the /p + 3.17 Up cut-off is, in a very good 
approximation, a single active electron-effect, meaning 
that at a given laser intensity only a single electron, 
namely the one that is next in the row for sequential 
ionization, contributes to the harmonic generation. Var- 
ious aspects of two- and many-electron effects on HOHG 
were studied in Refs. [1, l3,BH[l3, 111- However, to the 
best of our knowledge nonsequential double recombina- 
tion (NSDR) and the associated second plateau has not 
been revealed so far. 

The known cut-off at Ip -\- Z.VIUp can be explained 
within the so-called "simple man's theory" (see, e.g., the 
review Ref. [ist and references therein) : an electron is re- 
leased with vanishing initial velocity at the emission time 
te due to ionization by a laser field with a vector potential 
A{t) and an electric field E{t) — ~dtA{t). Thereafter, 
the electron moves freely in the laser field without be- 
ing affected anymore by the binding potential V{r), i.e., 
its momentum and position at time t > tc are given by 
p{t) = A{t)-A{t^) andr(t) = a{t)-a{Q-A{U){t-Q, 
respectively, where a.{t) — J dt' A{t') is the excursion of 
a free electron in the field (atomic units H,m, |e| , Aireo = 1 
are used unless noted otherwise) . In order for a harmonic 
photon being emitted, the electron has to revisit the ion 

at some time t^- > tc, that is, r{t^) — 0, since only then 
its overlap with the ground state — and thus the recom- 
bination probability — is appreciable. The energy of the 



emitted harmonic photon is given by fl = p^(ir)/2 + /p. 
Searching the pairs {tc, tr) for which r{t,-)^Q andp^{tr)/2 
is maximum, leads in the case of a constant-amplitude, 
linearly polarized laser field in dipole approximation, e.g., 
A{t) = AezSinivt, to max(p^(tr)/2) = 3.17C/p where 
Up = iV4. 

Let us now consider a two-electron atom or ion where 
the two electrons are freed by sequential ionization. The 
possible kinetic return energies — p'^{t^)/2 are shown 
in Fig. [H Each of the two electrons has a maximum 
return energy 3.17 Up. If both electrons moved along 
the same trajectory in the continuum and recombined 
together, emitting a single photon, one would expect a 



HOHG cut-off at 2-3.17C/o 
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with /A the ion- 



ization potential for the zth electron. However, electron 
repulsion renders this process extremely unlikely since 
the two electrons would have to "sit on top of each other" 
for about half a laser c ycle (which is even less likely than 
"collective tunneling" [I4I). Instead, if the electrons are 
emitted at certain times during subsequent or next but 
one half cycles, there are trajectories that do not cross 
before the recombination event, so that their mutual re- 
pulsion plays only a minor role (and is compensated by 
the attraction of the ionic potential). If, with respect 
to the first electron, the second electron is emitted dur- 
ing the subsequent half cycle, the sum of return energies 
around its maximum value is shown in Fig. [T] (indicated 
by "2 -I- 1"). For the constant-amplitude pulse one finds 
4.70 ?7p. If, however, the delay between the emission of 
first and second electron is greater than half a cycle the 
sum of the return energies can be even higher, as clearly 
visible from the values labelled by "3-1- 1" in Fig.[TJ The 
HOHG cut-off for NSDR is then expected at 



max(riNSDR) = 5.55 Up 
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In the following, we will show that this cut-off can be 
clearly identified in ab initio solutions of the two-electron, 
time-dependent Schrodinger equation (TDSE) and using 
an extended strong field approximation (SFA). 

First, we employ a widely used one-dimensional (ID) 
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FIG. 1: (color online). Return energies in units of Up vs the 
return time ti in laser cycles for a laser field A{t) = Ae^ sinuit. 
Electrons emitted during the first half cycle (1, black), the 
second (2, cyan), the third (3, red) etc. reach their maximum 
return energy 3.17 f/p during the subsequent half cycle. At 
later returns only lower values are achieved. For two electrons 
returning at the same time but emitted at different times, the 
values indicated by "2 + 1" and "3 + 1" are obtained (the sum 
"3+2" equals "2+1" in height and is omitted in the plot). 



model He-atom [IJ] to study NSDR on an ab initio TDSE 
level. The Hamiltonian in dipole approximation reads 
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with V{xi) 
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+ y(x,) )+Wixi-X2), 
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Ccij W{xi - X2) = [{xi - 

and pci the canonical momenta 
of the two electrons, i.e., Xi = pd + A{t). The soft- 
core parameters eci, £00 can be tuned in such a way 
that the model ionization potentials Ip'^ equal the real 
ones of the 3D He atom. The actual values were Cei = 
0.5 (which yields the correct /p^-* = 2.0 for He+) and 
£00 = 0.329 (so that = 0.904) The TDSE 

i9t^'(xi, a;2, i) = H{t)'^{xi,X2,t) was solved on a xi,X2- 
grid {Ax = 0.2) using a split-operator Crank-Nicolson 
(Peaceman-Rachford) approach (At — 0.075), starting 
at t = with the spatially symmetric spin-singlet ground 
state wave function. The n-cycle laser pulse was of the 
form A{t) = A sin^ {cut /2n) sin Lut ior < t < nT and 
zero otherwise (with T = 2it/lo). The harmonic spectra 
are calculated from the modulus-square of the Fourier- 
transformed acceleration a[t) [l^. 

Results for HOHG spectra obtained from the He model 
atom exposed to a n = 6-cycle laser pulse with w = 
0.0584 and A = 3.417 (/ = 1.4 • 10^^ Wcm'^) are shown 
in Fig. [21 Besides the result for the fully correlated sys- 
tem with both electrons active, the corresponding spectra 
for a frozen, inner electron (labelled "SAE" , drawn blue) 
and He^ (drawn red) are shown. For the single active 
electron (SAE)-calculation, a frozen Hartree-Fock poten- 
tial (leading to the proper ionization potential for the 
outer electron = 0.904) was used. Not surprisingly. 



FIG. 2: (color online). HOHG spectra for a n = 6-cycle laser 
pulse with uj = 0.0584 and A = 3.417 (/ = 1.4 • 10^^ Wcm"^). 
The spectrum for the He-model with both electrons active 
(drawn black) shows a second plateau. The corresponding 
single active electron (SAE)-result and the spectrum obtained 
from the He^-ion are also included. The arrows indicate the 
expected cut-off positions (see text). 



the SAE spectrum agrees well with the two active elec- 
tron (TAE)-result throughout the "usual" plateau (whose 
cut-off 3.17 Up + Ip'' pJI] is indicated by arrow 1) because 
at the chosen laser intensity it is the outer electron that 
contributes most to HOHG. As expected, the plateau of 
the He"'"-spectrum is several orders of magnitude lower 
in efficiency. The He"*" HOHG-cut-off is expected at 
3.17 f/p + 4^' (indicated by arrow 2). Neither the SAE- 
nor the He+-spectrum show a second plateau. The TAE- 
result, however, confirms our above considerations. Ar- 



rows 3 and 4 indicate the positions 4.70 C/p+/p 
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Pp' and 



/p , respectively 17]. Around arrow 3 a 



5.55C/p + 4^^ 

qualitative change occurs: for lower harmonic frequencies 
the spectrum displays a rich interference structure be- 
cause many "quantum trajectories" [isj contribute. The 
situation changes for the harmonic orders where only the 
classical solutions with t^/T G [3,3.5] survive. There are 
only two such solutions (and finally, at the cut-off, only 
one), resulting in a much less jagged spectrum between 
arrows 3 and 4. 

Switching-off electron-electron correlation in the two- 
electron TDSE-code is equivalent to the simulation of two 
independent He"'"-ions. For the latter, NSDR is impossi- 
ble, showing that electron correlation is clearly required 
for NSDR to happen since otherwise the two electrons 
cannot emit their total energy in a single photon. This 
is well known for the inverse process, i.e., single photon 
double ionization. 

The classical trajectories directly emerge in the TDSE 
quantum results if a time-frequency analysis is per- 
formed. To that end a window is applied to the spec- 
trum, selecting a certain frequency interval. The result is 
Fourier-transformed back, leading to a complex quantity 
a(yi, t) containing information about when the selected 
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FIG. 3: Time-frequency analysis of the HOHG spectrum in 
Fig. [2] showing the color-coded contour plot of logj^Q |a(f2, t) |^. 
The relevant classical solutions for the 6-cycle sin'^-pulse are 
superimposed. The possible classical recombination energies 
of a single electron Q = E^+Ip^^ are drawn white. The sums of 
the two highest classical return energies Eii , Er2 plus the two 
ionization potentials, Q. — X^i [-^ri + /p''], are superimposed in 
black for the recombination times tr of interest. 

harmonics are emitted. The time-frequency analysis of 
the TAE-spectrum of Fig. [5] is shown in Fig. [31 The 
relevant classical "simple man's" solutions are superim- 



posed (analogous to Fig. [T]). The recombination energies 
+ /p^'' of individual electrons are drawn white. In 
the frequency intervals of interest, the sum of the two 
classical highest return energies (plus the ionization po- 
tentials) are superimposed in black. The excellent agree- 
ment between "simple man's" and ah initio quantum re- 
sults shows that NSDR is indeed the mechanism behind 
the second plateau. 

The instants of ionization of the two electrons can be 
controlled by applying attosecond pulses of photon en- 
ergies close to the ionization potentials. In that way we 
were able to enhance the strength of the NSDR plateau 
in Fig. [2] by three orders of magnitude. More details will 
be given in a forthcoming publication. 

One may object that a ID model He overestimates 
the efficiency of the NSDR process because of a reduced 
wave-packet spreading while the electrons are in the con- 
tinuum, or other effects. The remainder of this Letter 
is therefore devoted to NSDR from an SFA perspective 
where the process can be studied in full dimensional- 
ity. Within the SFA the dipole in polarization direction 
^NSDRj-^) responsible for HOHG via NSDR can be derived 
similarly to the case of nonsequential double ionization 
[3 and reads cF.^^^it) = Y.l=i dT^^it) with 



t t' t' 
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The interpretation of (t) is as follows: electrons 

1 and 2 start both in their ground states l^i), |$2} 
with binding energies fi, £2, are dislodged by the elec- 
tric field E at times ti and t2, respectively, interact at 
time t' > ti, t2 via V12, and recombine at time t > t' . The 
Volkov time evolution operator uf^\ i — 1,2 governs the 
propagation of electron i in the laser field (without any 
Coulomb interaction). We evaluate ([T]) using hydrogenic 
ground states with energies £i = —0.9 and £2 — —2.0, 
the expansion of the lerigth gauge Volkov-propagators in 
Volkov- waves (see, e.g. [l3|), assuming a contact- type in- 
teraction V12 = 5{ri)5{ri — r2) and employing the 
standard saddle-point integration over intermediate elec- 
tron momenta [3[. In the case of a contact interaction, 
the triple time integral in the final expressions for the 
d^f^^{t) can be disentangled into a double time integral 
over the interaction time i', the ionization time i^, and 
an integral over the remaining ionization time tj^i that 
needs to be calculated only once for all t', substantially 
reducing the numerical effort. 

Figure 3] shows spectra obtained by Fourier- 



transforming (i^^^^ (i) for (a) the same laser parameters 
as in Fig.[2]and (b) a six-cycle pulse of twice the frequency 
and A = 2.0. A second plateau is obtained in both cases 
although with very different efficiencies. The SFA result 
for |d^^^^ (ri)pri^ is, owing to the wave packet spreading 
of the two electrons, proportional to uP' and uj^ in one and 
three dimensions, respectively. As a result, the level of 
the NSDR plateau in the low frequency case of Fig. |4li is 
much lower than in Fig. [2l The wave packet spreading is 
substantially reduced for the doubled frequency, leading 
to much more efficient NSDR in Fig. 0)3. The probabil- 
ity Pnsdr of NSDR to occur is Pnsdr = ^'hc^ho+-Poc 
where Prc, -Pho+ are the probabilities for HOHG from 
He, and He"^, respectively, and Poo is the probability for 
the emission of a single photon due to electron-correlation 
(instead of two lower-energy photons). Because in NSDR 
both electrons return to the nucleus at the same time, 
electron-electron interaction is very likely, i.e., Poc — 1 is 
expected, and confirmed by both the ID TDSE and 3D 
SFA results. The level of the NSDR plateau can thus be 
estimated by simply multiplying the levels of the SAE 
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FIG. 4: (color online). High-order harmonic spectra ob- 
tained by Fourier-transforming the SFA-NSDR expression for 
the dipole d^^^^'{t) (black). For comparison, the single ac- 
tive electron-result for He (blue) and He^ (red) are included; 
(a) same laser parameters as in Fig. [5] (b) lj = 2 ■ 0.0584, 
A = 2.0. Classically expected cut-offs are indicated by arrows: 
(1) |£-i| + 3.17f/p, (2) \£2\ + 3.17(/p, (3) \£i\ + \£2\ + 4.70[/p, 
(4) \£i\ + \£2\ + 5.55t/p, and (5) \£i\ + \£2\ + 2 ■ 3.17?7p. 

He and He+ plateaus. The unlikely events of simultane- 
ous tunneling ionization (ti = ^2) and motion of the two 
electrons along the same trajectory are not inhibited in 
the SFA dipole (P). As a consequence, the SFA-NSDR 
spectrum extends up to -I- |£2| + 2 • 3.17J7p, a fact 
which, however, is irrelevant for the determination of the 
relative strength of the second plateau. 

In conclusion, the existence of nonsequential double re- 
combination in intense laser fields was revealed. In this 
process, two electrons are freed by the laser, move in the 
continuum, and are driven back to the ion by the laser, 
where they recombine together upon emission of their 
joint energy plus their ionization potentials as a single 
photon. Nonsequential double recombination may thus 
be viewed as inverse single photon double ionization. It 
manifests itself in a second plateau in high-order har- 
monic generation spectra. This was demonstrated both 
by solving the time-dependent Schrodinger equation of 
a two-electron model atom numerically and using the 
strong field approximation. The efficiency of the pro- 
cess is approximately given by the product of the effi- 
ciencies for high-order harmonic generation in He and 
He"*". It may be substantially enhanced by aiding the 
desired ionization times with the help of xuv attosecond 
pulses. The main features of the harmonic spectra, in- 
cluding the position of the new cut-off, were explained 
in terms of classical trajectories. Finally, it should be 
noted that, instead of simultaneous recombination, the 
two returning electrons may undergo elastic or inelastic 
(2e, ne)-scattering processes such as, e.g., nonsequential 
ionization involving 2+n electrons or higher-order above- 



threshold ionization. 
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